KUNS-2211 



Instability of Small Lovelock Black Holes in Even-dimensions 

Tomohiro Takahashi and Jiro Soda 
Department of Physics, Kyoto University, Kyoto, 606-8501, Japan 

(Dated: October 4, 2009) 

Abstract 

We study the stability of static black holes in Lovelock theory which is a natural higher dimen- 
sional generalization of Einstein theory. We derive a master equation for tensor perturbations in 
general Lovelock theory. It turns out that the resultant equation is characterized by one functional 
which determines the background black hole solutions. Thus, the stability issue of static black 
holes under tensor perturbations in general dimensions is reduced to an algebraic problem. We 
show that small Lovelock black holes in even-dimensions are unstable. 
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I. INTRODUCTION 



It is well known that string theory can be formulated only in ten dimensions. Hence, it is 
necessary to reconcile this prediction with our real world by compactifying extra-dimensions 
or by considering braneworld. Intriguingly, in the context of the braneworld with large 
extra-dimensions, black holes could be created at the TeV scale Hence, the stability of 
higher dimensional black holes becomes important since these black holes could be produced 
at the LHC if the spacetime has larger than six dimensions. 

The stability of higher dimensional black holes has been an active topic since the seminal 
papers by Kodama and Ishibashi It is important to study various black holes in Einstein 
theory because black holes produced at the LHC are expected to be charged or rotating. 
A numerical study of charged black holes has been done |3|. To investigate the stability of 
rotating black holes, a group theoretical method is developed j4|. The method is used to 



study the stability of squashed black holes 



6j and 5-dimensional rotating black holes I? | 



The stability of rotating black holes in more than 5-dimensions is also studied 110|. 
It is also important to consider the stability of black holes in more general gravitational 
theories because black holes are produced at the Planck scale where Einstein theory would 
be no longer valid. In fact, it is known that Einstein theory is merely a low energy limit 



of string theory 



In string theory, there are higher curvature corrections in addition 



to Einstein-Hilbert term 



Thus, it is natural to extend gravitational theory into those 



with higher power of curvature in higher dimensions. It is Lovelock theory that belongs to 



such class of theories 121. In Lovelock theory, it is known that there exist static black hole 

n 

solutions [13|]. Hence, it is natural to suppose black holes produced at the LHC are of this 
type jl^. Thus, it is important to study the stability of these Lovelock black holes. 

In the case of second order Lovelock theory, the so-called Einstein-Gauss-Bonnet theory. 



the stability analysis under tensor perturbations has been performed [15| (see also an earlier 
work [l^). The analysis has been also extended to the scalar and vector perturbations jirj l. 
It is shown that there exists the scalar mode instability in five dimensions, the tensor mode 
instability in six dimensions, and no instability in other dimensions. In the case of third order 
Lovelock theory, the stability analysis of Lovelock black holes under tensor perturbations 
has been done by us 1^ . We have shown that there is the instability for small black holes in 
eight dimensions. Although third order Lovelock theory is the most general theory in seven 
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and eight dimensions, it is not so in more than eight dimensions. For example, when we 

consider ten dimensional black holes, we need to incorporate fourth order Lovelock terms. 

Indeed, when we consider black holes at the LHC, it is important to consider these higher 

□ 

order Lovelock terms [19|. Hence, in this paper, we study the stability of black holes in any 
order Lovelock theory, namely, in any dimensions. 

The organization of this paper is as follows. In Section [Tll we review Lovelock theory and 
explain a graphical method for constructing Lovelock black hole solutions. In Section IIIIl 
we present a master equation for tensor perturbations in the background of Lovelock black 
holes and reveal its universal structure. In Section HVl we examine the stability of Lovelock 
black holes with the method developed previously 18|. Finally, we summarize our results 
in Section IVl 

II. LOVELOCK BLACK HOLES 

In this section, we review Lovelock theory and introduce a graphical method to obtain 
asymptotically flat black hole solutions. 

In I2I, the most general symmetric, divergence free rank (1,1) tensor is constructed 
out of a metric and its first and second derivatives. The corresponding Lagrangian can be 
constructed from m-th order Lovelock terms 

r rAi(Ti--AmO-,„ p PlKl p PmK.m ( -\\ 

™ 2m Pl«l--PmKm^AiCTi ' ' ' -^^XmCr^n 1 \^ ) 

where Rxa^'^ is the Riemann tensor in D-dimensions and S^^^^'"^"^!^"^ is the generalized totally 
antisymmetric Kronecker delta. Then, Lovelock Lagrangian in D-dimensions is defined by 

k 

L = ^ ^ CffiCyji , (2) 
m.=0 

where we define the maximum order k = [{D — l)/2] and are arbitrary constants. Here, 
[z] represents the maximum integer satisfying [z] < z. Hereafter, we set cq = — 2A, ci = 1 
and Cm = O'm./'m- > 2) for convenience. Taking variation of the Lagrangian with respect 
to the metric, we can derive Lovelock equation 

k 

n _ r'l' _ XXf _ \ ^ XvXiai - XmO-m ID PlUl ID PniKm ( 0\ 

m=l 
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As is shown in 13j, there exist static exact solutions of Lovelock equation. Let us consider 
the following metric 



ds^ = -f(r)dt^ + -— + r'^'jijdx'dx^ 
f{r) 



(4) 



where is the metric oin = D — 2- dimensional constant curvature space with a curvature 
fi;=l,0 or -1. Using this metric ansatz, we can calculate Riemann tensor components as 

f 



R 



tr 



tr 



(5) 



Substituting ([5]) into (IHj) and defining a new variable ip{r) by 



/(r) = K, — r'^ip{r) 



(6) 



we obtain an algebraic equation 



m=2 



'2m-2 



m 



l[{n-p)\r 



p=i 



2A 



/i 



n(n 



(7) 



In (I7j), we used n = D — 2 and is a constant of integration which is related to the ADM 
mass as |20j: 



M 



(8) 



r((n+l)/2) ' 

where we used a unit IGnG = 1. 

From ([7]), it is easy to see that /(r) has many branches. In this paper, we want to 
concentrate on asymptotically fiat spherically symmetric, k, = 1, solutions with a positive 
ADM mass /i > because such black holes could be created at the LHC We also assume 
that Lovelock coefficients satisfy 



A = 0, a™ > 



(9) 



for simplicity. For example, consider n = 5 for which Eq.Q becomes third order algebraic 
equation. Though we have a formula for solutions of Eq.([7]) in this case, the roots are too 
complicated in general. Hence, we use a graphical method illustrated in FigHJ Because 
of the conditions the function is monotonic for positive ip. From ([7]), we see the root 
behaves ip ~ fi/r^~^^ or /(r) ~ 1 — fi/r"-"^ as r — oo. Thus, the asymptotically flat solutions 
belong to the branch where ijj is always positive. 



a3(n-l)(n-2)(n-3)(n-4)\|fV3 
^ +a2(n-l)(n-2)\|/V2+\|/ 


/ u/r"+l 


/ 1 
/ I 
/ 1 
/ I 
/ 1 
/ \ 








1 V(r) ^ix) ^ ^ 

i i 

1 asymtotically AdS branch asymtotically AdS branch 
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FIG. 1: We illustrate a graphical method for n = 5 case. In this case, third order Lovelock theory 
is a complete theory. Therefore, W\^\ in ([7]) reads a cubic polynomial. In this figure, three roots 
are depicted. Among these roots, only ^ > one corresponds to an asymptotically flat solution. 

Note that the horizon radius of the asymptotically flat solution is characterized by 
fi^n) = 0. From ([6]), we have a relation ipn = 'ip{f'H) = ^/^"h- Using this relation and 
([7]), we obtain an algebraic equation 

W[^h] = f^^^H^'^^' . (10) 

This determines ipn and hence vh- In FigjJ], we present a graphical method to solve Eq. (fTO!) . 
It is obvious from Figl2]that oo > r > th corresponds to < < ipH when /(r) describes 
an asymptotically flat solution. It is also obvious that ipH becomes larger as /i is smaller. 
Note that there may be no horizon if is too small (for example, in seven dimensions, there 
is no horizon if < /i < 803 — Oog). 

Using the metric ansatz (jl]), the Kretschmann scalar Rabcd.R"'^'^'^ is 

RabcdR""""' = f" + 2n^ + 2n{n - l) ^"" ~ P\ 

Then, this solution has curvature singularities at r = or where derivatives of /(r) diverges. 
For asymptotically fiat spherical symmetric solutions, the conditions |1 — /| < 00, |/'| < 00 
and I/" I < 00 are satisfied except for r = 0. Therefore, there is a curvature singularity only 
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FIG. 2: For n = 5 case, this figure explains a method for calculating ipn or rn graphically. 



at r = in the cases we are considering. And asymptotically flat solutions have a horizon if 
the parameter /i is sufficiently large. So, this solution doesn't have a naked singularity and 
describes a black hole with a mass M defined in ([8]) and no angular momentum. 

III. MASTER EQUATION FOR TENSOR PERTURBATIONS 

To analyze the stability, it is convenient to decompose the metric under the symmetry of 
n-dimensional symmetric space. There are scalar, vector, and tensor modes. In this paper, 
we study tensor perturbations. In this section, we do not restrict parameters. 

We consider tensor perturbations around the solution (jlj) 

6gab = , Sgai = , 5gij = r'^(l){t,r)hij{x') , (11) 

where a,b = {t, r) and r) represents the dynamical degrees of freedom. Here, hij are 
defined by 

V'^Vfe/i.j = -ihj , = , f% = 0. (12) 

Here, V* denotes a covariant derivative with respect to 7jj. Here, the eigenvalue is given by 
7 = + n — 1) + 2, (^ = 2, 3, 4 ■ ■ ■ ) for k = 1 and negative real number for k = —1, 0. 
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Perturbing the Lovelock equation ([3]), we can get the first order perturbation equation 



= = - 



k 



ar. 



/ V 2(m+l) MPi'^i'"P™''"i 

771=1 



'Am-icr 771-1 



Pm — 1 '^m — 1 



5i? 



' A 771 CT 777 



pm 



(13) 



where 5Rab^'^ represents the first order variation of the Riemann tensor. Since we are con- 
sidering tensor perturbations, apparently we have 5Qa = 5Qa = 0. To calculate 5Ql, we 
only need SRu^^ , 5Rri'^ and dRij^K The reason is simple. We consider, for example, the 
term proportional to SRtk"^ in 6Qi^ . This means \m = t or am = t in (fT3!) : and Pm 7^ t and 
Km 7^ t. Since ^pj^^.'.'.' consists of Kronecker delta, there must exist i {1 < i < m — 1) such 
that Pi = t OT K£ = t. Then or ag must be t because of background Riemann tensor (Q. 
So, for example. Am and \e become t. Noticing that Spl^l''' is totally antisymmetric, we see 
Splnl--- = when Am = A^. Thus, the coefficient of SRtk "^ in 6Qi^ must be 0. For the same 
reason, coefficients of other than 5RuK SRri^^ and 5Rij^^ vanish. From the metric ansatz 



the necessary components can be deduced as 



SRu'' = 

5RrP = 
'J I O-Tljm — 

Thus, using the following relation 



2/ 



4 r 



n — 2 f , 2k — 7 



2 r 



2r2 



(14) 



"ikik2-k^-ik^"l^ — "'')"ikik2- km-i ' 



m I'm 



(15) 



we can calculate SQ-! as follows: 
/i(r 



rn—2 



-,11—2 



{6Ru'' + 5RrP) + 



h 



where 



h{r) = r"-2 + 



m=2 



n - 2)r"-3 

f'h , fh , fh'\ ^, , (2k-7)/i' 

H ; ^ — ' 

2 r 2 y ^ 2(n - 2)r 



^2m-2 

IT-' ( n - 

p=2 



(16) 



X <^ -(m - l)r/' + {n- (2m - !))(«: - /) 



(17) 
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It is useful to note that the universal function h is related to W as 

"'^-1 dWH 



h{r) 



d 
dr 



n 



1 di) 



(18) 



Separating the variables 4>{r,t) = x{r)e'^^, we can deduce the master equation from (fT3|) 
and f|T6|) as follows: 



-fx 



h r J [n — 2)r n 



(19) 



Here, we should notice that we have not used information about Lovelock coefficients and 
/(r). Hence, this master equation itself is very general, although we will consider specific 
situations in the next section. Moreover, this equation is the same as the master equation 
in |l^ except for the concrete form of h{r). Hence, the method used in 18| can also be 
applicable. 



IV. EVEN-DIMENSIONAL SMALL BLACK HOLES ARE UNSTABLE 

In this section, we first clarify condition for the stability of black holes. Then, we demon- 
strate that small black holes are always unstable in even-dimensions. 



A. Condition for Stability 

In this section, we assume the conditions (Q are satisfied. Then, there exists an asymp- 
totically flat branch. And we also assume spherical symmetry and positivity of the mass, 
i.e., K = 1 and /i > 0. 

As we will soon see, the master equation ( fT9l) can be transformed into the Schrodinger 
form. To do this, we have to impose the condition 

h{r) > , (for r > riy) . (20) 

In fact, this is necessary for the linear analysis to be applicable. In the case that there exists 
tq such that h{ro) = and tq > th, we encounter a singularity. In fact, using approximations 
h{r) ~ h'{rQ){r — Vq) = h'{rQ)y, f{r) = /(ro) and r = r^, we can show the solution near Vq 
is given by x ~ Ci + C2logy, where Ci and C2 are constants of integration. The solution is 
singular at y = for generic perturbations. The similar situation occurs even in cosmology 
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with higher derivative terms |2l|,|22|. In those cases, this kind of singularity alludes to ghosts. 
Indeed, there is a region h{r) < outside horizon where kinetic term of perturbations has 
wrong sign. Hence, we simply assume these black hole solutions are not allowed. 



When the condition fl2U]) is fulfilled, introducing a new variable \E'(r) = x{^)'^\/h{r) and 
switching to the coordinate r*, defined by dr* /dr = 1//, we can rewrite Eq.( |T9l) as 



dr*'^ 

where 



+ \/(r(r*))^ = -tu^^ = , (21) 



, (2hi — ■y)f dlnh 1 ^d / d r-\ 

is an effective potential. 

For discussing the stability, the " S- deformation" approach is useful {2!, Q . Let us define 
the operator 

acting on smooth functions defined on / = (r|^, 00). Then, (1211) is the eigenequation and E 
is eigenvalue of Ti. We also define the inner products as 



((^1,(^2) = j^V*iV2dr* . (24) 
In this case, for any ip, we can find a smooth function S such that 

{ip,nip) = Ji\Dip\' + V\ip\')dr*, (25) 

where we have defined 

D = :f-^+S, V = V + f^-S\ (26) 
dr* dr 



Following [15[, we choose S to be 



S = -f^\n{r^/h) . (27) 



Then, we obtain the formula 

roo \,„\2 



n^) ^ I |«.r^*. + (2. - -,) £ J^^* . (28) 
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Here, the point is that the second term in (128!) includes a factor 2^: — 7 > 0, but h does not 
include 7. Hence, by taking a sufficiently large 2k — 7, we can always make the second term 
dominant. 

Now, let us show that the sign of dlnh/dr determines the stability. If d\nh/dr > on 
/, the solution (jl]) is stable. This can be understood as follows. Note that 2k — 7 > 0, 
then we have > for this case. That means {ip^Tiip) > for arbitrary if ii dlnh/dr > 
on /. We choose, for example, if as the lowest eigenstate, then we can conclude that the 
lowest eigenvalue Eq is positive. Thus, we proved the stability. The other way around, if 
dlnh/dr < at some point in /, the solution is unstable. To prove this, the inequality 

^ > E„ (29) 

is useful. This inequality is correct for arbitrary if. If dlnh/dr < at some point in /, we 
can find if such that 

f°° lifl^ dlnh , 

/ . ^ ^dr < . (30) 

(n-2)r dr ^ ^ 

In this case, (1281) is negative for sufficiently large 2k, — 7. Then, the inequality (129!) implies 
Eq < and the solution has unstable modes. Thus, we can conclude that the solution is 
stable if and only if dlnh/dr > on I. 

To summarize, we need to check the sign of dh/dr outside the horizon r > rn to inves- 
tigate the stability. Note that the condition h{r) > is necessary to have healthy black 
holes. 



B. Instability of Small Black Holes 

Here, we assume h{r) is positive outside the horizon. Otherwise, such black holes are 
already sick. Hence, what we have to do is to check the sign of dh/dr. 

In order to see the sign of dh/dr, it is very useful to express it as function of tp instead 
of that of r. Using ([71) and its derivative, we obtain 

{d^Wm ^' = -{n + l)-^^ = -{n + 1)^ . (31) 

The above formula can be used to eliminate ip' in (1T8|) . The result reads 

M.) ^ (a,w - '1^"^) . (32) 



n-1 d^W 
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Similarly, dh/dr can be written as 



dr 



+ 



n - 1 d^W 
n + 1)2 {d^WdlW - {dlWf] 



(33) 



n - 1 {d^Wf 

Since is a polynomial function of ip and d^W is positive, we can determine the sign 

of dh/dr by examining the sign of polynomial determined by the above formula. Thus, the 
stability problem has been reduced to an algebraic one. 

Substituting the explicit form of W[ip] into this equation, we obtain 

^ = r»-= 3 ^ (34) 

*■ (i + EL.Kjn^ri^-p)}'/-]) 

Here, the lowest and leading term of L{x) is 
L{x) = (n - 2) + ■ ■ ■ 

4 f2k-2 ^ 

+ - 1)=^ I n - J - - - - l))a;^*^-^ . (35) 

We note that the highest order k = [{D — l)/2] is related to the dimensions as n = 2A; — 1 
in odd- dimensions and n = 2k in even-dimensions. In odd-dimensions, the leading term 
disappears. Hence, we cannot say anything in general. Hence, we consider only even- 
dimensions. 

Let us examine the sign of L{x) (x > 0). If n = 2k, the coefficient of the lowest term 
is positive and that of the leading one of fl35l) becomes negative. Therefore L{x) > near 
X = and L{x) < for large x. This means that there exists roots of L{x) = because L{x) 
is continuity function; we define Xq as the lowest positive root. If tpn < Xq, then L[ip] > 
for < tp < ipH, and hence we see dh/dr > for r > ru- While, if ipn > xq, then there 
exists a region L[ip] < in xq < ip < ipn, and thus there exists an area dh/dr < in r > r^j- 
Therefore, black holes are stable if ipn < xq and unstable if tpn > xq. Considering the result 
that ipH becomes larger as /i, which is related to ADM mass, becomes smaller, we conclude 
that there exist critical mass below which black holes become unstable. 

Summarizing this section, there is critical mass below which black holes are unstable 
in n = 2k. Thus, we have shown the instability of small Lovelock black holes in even- 
dimensions. 
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V. CONCLUSION 



We have studied the stabihty of static black holes in Lovelock theory which is a natural 
higher dimensional generalization of Einstein theory. We derived a master equation for 
tensor perturbations in general Lovelock theory. It turned out that the resultant equation is 
characterized by one functional W[ip] which determines the background black hole solutions 
through Eq.©. The stability issue of static black holes under tensor perturbations in general 
dimensions has been reduced to an algebraic problem. We have shown that there exists the 
instability of Lovelock black holes with small masses in even-dimensions. Remarkably, the 
instability is strong on short distance scales. Curiously, the similar instability also appears 



in the Gauss-Bonnet cosmology 



23|. 



ack 



17|, if 



We have done the analysis of tensor perturbations and found the instability of smal 
holes in even-dimensions. Taking into account the results in Gauss-Bonnet theory 
we extend the analysis to scalar and vector perturbations, we can expect that small black 
holes are unstable also in odd- dimensions. If so, we need more profound understanding of 
the instability in Lovelock theory. 

It is also interesting to investigate the fate of the instability. As the instability is stronger 
for higher multipole orders i, the resultant geometry would be weird. This issue is very 
important because black holes lose their mass due to the Hawking radiation and eventually 
become unstable. 

Related to the above, we have to find meaning of the universal function h{r). As was 
shown in this paper, this function governs the stability of black holes. Therefore, if h(r) 



has, for example, thermodynamical meaning, the relation between thermodynamical 
and dynamical instability might be revealed. 
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